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Abstract 

In this chapter, we show how the MEP hypothesis may be used to 
build simple climate models without representing explicitly the energy 
transport by the atmosphere. The purpose is twofold. First, we assess the 
performance of the MEP hypothesis by comparing a simple model with 
minimal input data to a complex, state-of-the-art General Circulation 
Model. Next, we show how to improve the realism of MEP climate models 
by including climate feedbacks, focusing on the case of the water- vapour 
feedback. We also discuss the dependence of the entropy production rate 
and predicted surface temperature on the resolution of the model. 

1 Introduction 

Although it is not straightforward to define what climate is precisely, 
one may suggest that what we call the climate system is made up of the 
atmosphere, the oceans, the cryosphere, the biosphere and the lithosphere 
[39] . The different components interact in various ways, and their relative 
importance depends on the question asked. For instance in numerical 
weather prediction, taking place on a timescale of a few days, the main 
dynamical component is the atmosphere and all the other components 
may be regarded as prescribed. On the contrary, the evolution of climate 
on very long timescales (of the order of tenths or hundreds million year) is 
essentially determined by the exchanges of carbon between the land, the 
oceans and the atmosphere. 

The distribution of surface temperature is of primary interest. It de- 
pends on a large number of factors, such as the composition of the at- 
mosphere (upon which the radiative energy exchanges depend), the cir- 
culation of the atmosphere and oceans, the ocean salinity, the presence of 
ice-sheets, the type of terrestrial vegetation cover,... State-of-the-art cli- 
mate models, usually referred to as General Circulation Models (GCMs), 
now include many of the above factors (the term Earth System Models is 
starting to emerge). 

However, not all this complexity is necessary to obtain a rough esti- 
mate of the temperature of a planetary atmosphere: perhaps the simplest 
approach is to balance the incoming solar radiation with the outgoing 
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planetary radiation. Again this can be done at various levels of accuracy, 
depending on the knowledge we have of the concentration of the radia- 
tively active constituents of the atmosphere (e.g. water-vapour and car- 
bon dioxide). Imposing a local radiative equilibrium is in fact misleading: 
latitudinal and vertical differential heating trigger atmospheric motions, 
which carry heat to mitigate the temperature gradients that would ex- 
ist at radiative equilibrium. The resulting energy transport term can be 
parametrized (for instance as a diffusion process with empirical diffusiv- 
ity) as a function of the temperature distribution, so that we can solve 
the model without resolving explicitly the motions of the atmosphere. 
Such models, consisting of a radiative model and a parameterization of 
the energy transport by the atmosphere are called Energy Balance Models 
(EBMs). Alternatively, one may solve the fluid dynamics problem and 
compute explicitly the velocity field: this is what GCMs do. The hierar- 
chy of climate models, ranging from simple EBMs to complex GCMs, also 
comprises the so-called intermediate complexity models (EMIC), which of- 
fer a variety of simplified representations of the atmospheric and oceanic 
circulation and other phenomena [33]. The main interest of EMICs is 
their relatively low computational cost, compared to GCMs, which make 
them particularly suitable for the study of palaeoclimates. Indeed, the 
timescales involved in such problems reduce the role of GCMs to sim- 
ulating snapshots. Both GCMs and EMICs require a certain amount 
of parameter tuning. This is sometimes a problem when studying past 
climates for which little data is available on which to base adjustment 
procedures, and even more so for other planetary climates, where many 
features differ tremendously from the terrestrial conditions on which the 
empirical parameterizations were tested. 

Nevertheless, the laws of physics remain the same when going back 
into time or out into the cosmos. The three branches of physics which 
play a fundamental part in setting the climate of a planet are radiation 
physics [IT], fluid dynamics |17[|38] and thermodynamics fil. One fun- 
damental principle which is always present, even in simple models like 
EBMs, is the first law of thermodynamics, because it describes the ex- 
changes of energy in a system. To energy exchanges are associated equi- 
librium temperature distributions. On the other hand, even in the most 
sophisticated climate models to date, the second law of thermodynamics, 
which also describes the exchanges of energy in a system but in a qual- 
itative rather than quantitative way, is not taken into account. When 
subgrid-scale parameterizations are involved, classical models may even 
violate the second law of thermodynamics 16 . It has also been suggested 
that spurious sources of entropy production could lead to a global cold 
bias in climate models [20]. Henceforth, a number of diagnostic tools 
emerged to study the thermodynamic properties of climate models [2 , 29 . 
Besides, postulating that the system chooses the steady-state with max- 
imum entropy production given certain constraints leads to a variational 
problem which has proved very efficient for predictive use. This is the 
so-called Maximum Entropy Production principle [24|[3l][34]. We shall 
not discuss here the theoretical foundations (or lack thereof) of this hy- 
pothesis (see 3,5,6, 13 ), but only its consequences for climate modelling. 
Hitherto, mainly two approaches have been developed. One point of view 
is that the MEP principle can be useful to select the value of adjustable 
parameters in empirical parameterizations from existing models, in an ob- 



jective way 19,22 23 25 37 . In the second approach, the purpose is to 



build simple climate models based on the MEP hypothesis for describing 



2 



unresolved processes. We shall present the latter approach in this chap- 
ter. After briefly reviewing earlier attempts (Sect. [5| we build a MEP 
climate model devoid of ad hoc assumptions and we show how to include 
feedbacks like the water- vapour feedback (Sect. [3]). The model is then 
tested for pre-industrial and Last Glacial Maximum conditions (Sect. [4]). 



2 The Paltridge model 



A typical one-dimensional EBM consists of a certain number of boxes, rep- 
resenting latitudinal zones, characterized by a single temperature. Each 
box receives energy from the outside in the form of solar radiation, and 
radiates back to space in the longwave domain. The difference of these 
two terms, which is usually called the radiative budget of the box, does not 
necessarily vanish: there are also energy exchanges with the neighbouring 
boxes due to atmospheric (and oceanic) transport of heat. Hence, for box 
i, the total energy budget reads 

where c Pi ,Ti,Ri and 7* denotes respectively the heat capacity, tempera- 
ture, radiative budget and atmospheric (or oceanic) convergence for box i. 
A radiative scheme provides Ri as a function of Ti'. e.g. Ri — £iS — eioTf 
where S is the solar constant, & represents the projection of the sur- 
face of the latitude belt onto the sphere centered on the sun, a is the 
Stefan-Boltzmann constant and e% the emissivity of the surface. In such 
a radiative scheme, the greenhouse effect is not taken into account. In 
contrast, there is no simple expression for 7^ which can be justified from 
first principles. A standard parameterization in this context is to assume 
a diffusion-like term, but there is no justification for this hypothesis and 
the diffusion coefficient has to be chosen empirically. 

Paltridge 35 suggested a model, with a more elaborate radiative 
scheme - involving in particular a cloud cover variable Oi in each box - 
than our above example, in which 7^ is not empirically parameterized as a 
function of the temperatures T , but instead satisfies a maximum entropy 
production principle. He postulates that the steady-state temperature dis- 
tribution Ti is such that the material entropy production rate a = ^\ ^f- 
is maximum, subject to the global steady-state constraint 7* = 0. At 
steady state, ji = —Ri and a is a function of the temperatures Ti only. At 
steady- state, the distributions of temperature, cloud cover, atmospheric 
and oceanic meridional fluxes obtained are in striking accordance with ob- 
servations. In spite of this apparent success, some major criticism remain. 
First of all, the planetary rotation rate is believed to be a major driver 
of the latitudinal distribution of temperatures, but it does not appear at 
all in Paltridge's model. Besides, it is clear that the principle does not 
hold in the case of a planet without atmosphere (see [9]). One may thus 
wonder if it is not pure coincidence that it seems to apply to the Earth's 
atmosphere 142]. Last but not least, there is no theoretical justification 
for the principle of maximum entropy production. 

The thread was taken up in a series of papers 10 12 36,46), verify- 



ing Paltridge's results in different variants of the original model, but the 
fundamental objections mentioned above remained unanswered. More re- 
cently, Lorenz [28] added some support to the idea that the agreement 
between the model and observations is not a coincidence, by showing that 
it gives acceptable results for Titan and Mars as well. The question of 



3 





1 


sa 






1 




T as 

T a A 




H 


±g l_l 







Figure 1: One grid cell of a two-layer MEP model. The surface layer has 
temperature T g and exchanges heat q (thick solid red arrow) with the overlying 
atmospheric layer of temperature T a . Both layers absorb solar radiation (thin 
solid yellow arrows) and emit and absorb longwave radiation (thin dashed blue 
arrows). The atmospheric layer exchanges energy with the surrounding cells: 
the convergence of the atmospheric heat flux is £ (thick solid red arrow). 

the independence with respect to the planetary rotation rate was also 
adressed by Jupp [2l] in a MEP model with a simple parameterization of 
atmospherics dynamics. Nevertheless, one fundamental concern remains: 
the Paltridge model and its variants still contain a large number of pa- 
rameterizations, ad hoc hypothesis and empirical coefficients, for instance 
in the radiative scheme, in the cloud parameterization or in the treatment 
of surface heat fluxes (maximum convective hypothesis). Is it possible to 
get rid of these potential biases to assess the intrinsic value of the MEP 
conjecture in the climate modelling framework ? This is the question we 
address in the next section. 

3 A simple MEP model with water- vapour 
feedback 

3.1 NEF Radiative scheme 

A possible strategy to assess the degree of coincidence in Paltridge 's re- 
sults may be to build a MEP model devoid of any ad-hoc parameter and 
assumptions. To that end, we suggest a new radiative scheme based on 
the Net Exchange Formulation (NEF), which only involves physical quan- 
tities (values of which are known a priori). Following [15] ? we introduce a 
two dimensional model with two layers: for each grid point characterized 
by a latitude and a longitude, there is a surface layer with a temperature 
T g and an atmospheric layer with a temperature T a . 

Each layer absorbs an amount of solar radiation (^g^ for the surface 
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layer and ^ as for the atmosphere) given by: 

V™ = (s(a g ) - s)(l - a g )£S, (2) 
= (s + a g s*)ZS, (3) 

where S is the solar constant, £ the projection of the cell area onto the 
sphere, a g the surface albedo, and the coefficients s, s* and s are adapted 
from the classical Lacis and Hansen scheme [26] : 

.-(«,) = 0.353 + ^-RA0 : A oz (Muq 3 )^ (4) 

1 — R*a g 

s A wv (Mu), (5) 

s* = ^((m+^|-^(Mm). (6) 

Here uo 3 , tfc represent respectively the vertically integrated ozone and wa- 
ter vapour density (including pressure scaling [44]) , M accounts for the 
slant path of solar rays, R r (Q and R* account for Rayleigh scattering 
in the atmosphere, and A oz ,A wv are absorption functions for ozone and 
water vapour. See 15,26,44 for details. 

The long-wave radiative exchanges can be written in a simple form 
using the Net Exchange Formulation 8 . The surface layer and the atmo- 
sphere exchange a net amount of energy through infrared radiation, 
while the surface and the atmosphere radiate respectively ^fip" and 
to space (see 15 for a derivation): 

= t(T g )oT2-t(T a )oTi, (7) 
¥™ = t(T a )alt, (8) 

where \i is the Elsasser factor arising from the angular integration, and 
t(T) = fi — J + °° Tjydi^j represents the emissivity of the atmo- 

sphere (B u is the Planck function). The transmission function r u de- 
pends on the vertical profiles of absorbing gases, pressure and tempera- 
ture: t v — exp ^— k iy (z)dz^j, where k u is the absorption coefficient, 

and H the total height of the atmosphere. To sum up, the only parameters 
required by the radiative scheme are the vertically integrated concentra- 
tions of water vapour u, carbon dioxide uco 2 (they determine k„), ozone 
uo 3 and the surface albedo a g . 

The steady- state condition for each box reads, for every grid point: 

*r+*fr-<-*L fl +c = o, (io) 
*ir = o, (ii) 

where £ is the horizontal convergence of atmospheric heat fluxes and q the 
surface to atmosphere heat flux. The total material entropy production 
is given by 

Nl at N lon , . 

*M{{T a , i3 ,T g , i3 }) -EE + ^ A - W 

\J-a,ij 1 g,i3 J-a^j J 

where Aij is the area of the grid cell in position and qij,Cij are 



functions of T a ,ij, T g ^j given by ( 10 )-( 11 ) . We are interested in the fields 
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that maximize <jm while satisfying the global constraint J]\ . AijQj = 0, 
which can be translated into an unconstrained variational principle using 
Lagrange multipliers. 

3.2 Different versions of the model 

The MEP model described in the previous section requires only physi- 
cal parameters as an input. In a first step, we compute the horizontal 
distribution of u (vertically integrated water vapour density) and uo 3 by 
linear interpolation of standard atmospheric profiles [32] (depending only 
on the latitude). To compare with the results of Paltridge, we also as- 
sume that the coefficients t(T) in Eqs. [7] are fixed, with a prescribed 
reference temperature T re f (dependent on the latitude) also computed 
from the standard profiles 32 (version vO in table [I]). However, the as- 
sumption of constant t(T) coefficient is very unrealistic: the shift in the 
Planck spectrum associated with a variation in temperature of the surface 
or atmospheric layer has a strong impact on the optical properties of the 
atmosphere. In version vl, we retain the dependence of the emissivity of 
the atmosphere on surface and atmospheric temperatures. Besides, fixing 
the profiles of water- vapour and ozone is also a restrictive hypothesis, es- 
pecially in view of potential applications to different climates for which 
standard profiles are not well known. As far as ozone is concerned, we can 
simply examine a version of the model in which we completely ignore ozone 
(version v2). For water- vapour, the situation is slightly more complicated: 
the atmospheric temperature is linked via the Clausius-Clapeyron relation 
to the water vapour content, which itself feeds back onto the temperature 
via the greenhouse effect. Yet, in the previous versions (v0-v2) of the 
MEP model, we kept fixed the absolute amount of water vapour in the 
atmosphere, independently of the temperature. In version v3, we fix the 
relative humidity RH = Ph 2 o / P S at(T). The vertically integrated den- 
sity of water vapour is related to the relative humidity, temperature and 
pressure profiles through: 

1 Mh 2 o f Ps uu „ d frr\ d P f1Q\ 
% 2 o = --w — / RH x P sat (T) — , (13) 

y WJ-air Jq P 

where Mh 2 o, M a % r are the molar masses of water and air, g is the gravity 
and P s the surface pressure. In our model with one atmospheric layer, we 
may assume that the relative humidity is uniformly distributed in each 
atmospheric cell, with a vertical extent equal to the scale height for water 
vapour. Relation (|13|) then becomes u* H<2 o ~ Mh 2 o /(gM a ir) x RH x 
P sa t(T) (version 3). The different versions are summarized in Table [I] 
The purpose of comparing these different versions of the model is at the 
same time to test the impact of reducing the quantity of input parameters 
(no T re f, no uo 3 ) and to improve the realism (Planck spectrum, water- 
vapour feedback). 

3.3 Water-vapour feedback and multiple steady 
states 

The physical quantities involved in the climate system are related in many 
ways, so that a change in one of these quantities can have an influence 
on another one, feeding back onto the original quantity, either moving it 
closer (negative feedback) or farther (positive feedback) from its initial 
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Model Version 


u 


Uo 3 


uco 2 




t(T) 


(Tpi)CC) 


(Tlgm — Tpi) 


MEP vO 


MC 


MC 


280ppmv 


T = 


Tref (MC) 


22.9 


-1.98 


MEP vl 


MC 


MC 


280ppmv 


T 


— TT 

— J-at J-g 


22.3 


-1.84 


MEP v2 


MC 





280ppmv 


T 


— TT 

— J-at J-g 


22.5 


-1.84 


MEP v3 


u*(T a ) 





280ppmv 


T 


— TT 

— J-a-, J-g 


19.9 


-2.9 


IPSL 






280 ppmv 






15.7 


-2.53 



Table 1: Different versions of the MEP model and the resulting global mean 
surface temperature for pre-industrial (PI) and last glacial maximum (LGM) 
climates, compared to GCM runs with the IPSL.CM4 model. "MC" stands 
for the integrated standard McClatchey profiles, and the angular brackets mean 
global average. See Sect. |3.2| for the definition of the different versions and Sect. 
[4] for the discussion of the results. 

value. A classical example of positive feedback is the water- vapour feed- 
back. If the temperature increases locally, the water vapour saturation 
pressure will increase so that more water (if available) may evaporate in 
the atmosphere, leading to stronger greenhouse effect and thus further 
increase of the temperature. Feedbacks of this sort can lead to multiple 
equilibria, bifurcations and hysteresis phenomena. For a given relative 
humidity distribution, equilibrium states with radically different temper- 
atures are simultaneously possible [41] . The water- vapour feedback has 
been shown to play a major part in important climate problems ES], ex- 
actly like feedbacks of different natures [27[|43| . Hence, it is essential to 
be able to represent them correctly in a climate model. In the context of 
MEP models, it was shown in LL4] that the ice-albedo feedback gives rise 
to multiple local maxima in the entropy production rate , corresponding 
to the multiple equilibria that appear in a traditional EBM (see also [2]). 
Here, we observe multiple local maxima of the entropy production rate in 
a certain range of solar constant and relative humidity. One great advan- 
tage of MEP is the small computational cost of maximizing a function as 
compared to integrating a complex differential equation. Of course this 
is no longer true if the function, or the submanifold on which to search 
for the maximum, becomes too complicated. Already, in the presence of 
multiple maxima, this difficulty has to be dealt with as the steady-state 
selected by the maximization algorithm may depend on the initial value. 
To avoid being trapped in an irrelevant state, several methods may be 
investigated. First it is possible to further restrict the manifold defined 
by the constraints to ensure that it contains only one local maximum 
of the entropy production. In the case of the water-vapour feedback in 
our two-layer model, solving the radiative balance for the whole column 
in terms of the atmospheric temperature may lead to several solutions. 
Selecting systematically one of them before computing the entropy forces 
the system to remain on the portion of interest in phase space. This is the 
technique that we use here. Alternatively, introducing the time dimension 
and assuming that at each time step, the system maximizes instantaneous 
entropy production with an additional term corresponding to time deriva- 
tives, it was suggested in 114 to use relaxation equations as a numerical 
algorithm to compute the final state (see also 45 ) . 
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Figure 2: Surface albedo a g in the IPSL model, for pre-industrial (left) and 
Last Glacial Maximum (right) conditions. 

4 Results: Present and Last Glacial Max- 
imum climates 

We compared the surface temperature distribution obtained from MEP 
with that obtained from a state-of-the-art GCM, the IPSL_CM4 model. 
The IPSL model is a coupled atmosphere-ocean model 30 used for the 
Fourth Assessment Report (AR4) of the Intergovernmental Panel on Cli- 
mate Change (IPCC) [18] . For pre-industrial climate, the forcings in the 
IPSL model are: pre-industrial greenhouse gas concentration (CO2=280 
ppm, CH4=760 ppb, N2O=270 ppb), insolation, coastlines, topography 
and land-ice extent. The surface albedo is computed from the IPSL_CM4 
pre-industrial simulation and used as a forcing for the MEP model (Fig. 
|] left). 

The surface temperature distribution obtained with the MEP model 
is represented in Fig. [3] along with the difference between the MEP model 
and the IPSL model. The global mean surface temperature for the MEP 
model is (Tn) = 22.9 ° C. By comparison, (Tpj) in the IPSL simulation 
is approximately 7°C lower (Table [T]); as Fig. [3] reveals, the major part 
of this difference comes from areas where the cloud cover is important, 
or elevated areas like the Antarctica. It is shown in 15 that a crude 
estimation of the effect of clouds and elevation suffices to explain the 
major part of the difference with the IPSL model. Figure [4] shows the 
meridional energy transport as a function of latitude for both the MEP 
model and the IPSL model for pre-industrial conditions. The agreement 
is remarkable given the simplicity of the MEP model. 

One advantage of the reformulation of the Paltridge model presented 
here is that due to the absence of ad-hoc parameters, it is possible to 
test the model on climates other than the Pre-industrial period. For 
instance, it is possible to change the surface albedo to take into account 
the variations of ice or vegetation extent. A time period which is largely 
documented and for which simulations with GCMs are available is the 
Last Glacial Maximum (LGM). It corresponds to the time during the last 
glacial period when the ice-sheets extent was maximum, roughly 21 000 
years ago El. At that time, large ice-sheets covered North America and 
Northern Europe, and the global mean temperature was approximately 

5 ° C lower than present. In the MEP model, it is only possible to take into 
account the effect in surface albedo due to the presence of the ice-sheets at 
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Figure 3: Left: surface temperature T g for pre-industrial conditions obtained 
with the MEP model (version vO). Right: Difference between the surface tem- 
perature T g in the MEP model and the IPSL model for pre-industrial conditions. 
Contour lines interval is 10 ° C, positive contours are drawn in solid lines, neg- 
ative contours in dashed lines and the null contour as a dotted line. 




Figure 4: Meridional energy transport as a function of the latitude in the MEP 
model (version vO, solid red line) and the IPSL model (dashed black line), for 
pre-industrial conditions. 
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Figure 5: Surface temperature difference between the Last Glacial Maximum 
and the Pre-Industrial, in the MEP model (left, version vO) and in the IPSL 
model (right). Contour lines space is 10 ° C, positive contours are drawn in solid 
lines, negative contours in dashed lines and the null contour as a dotted line. 



the LGM (Fig. [2] right), and not, for instance the associated topography 
effect. To ensure the comparison with the IPSL model is as direct as 
possible, we use a simulation where only the albedo effect is taken into 
account in the GCM. The resulting surface temperature difference between 
the LGM and the PI is shown in Fig. [5] for both models. The global mean 
difference is ~ — 2 C in the case of the MEP model and ~ — 2.5 ° C for 
the IPSL model. However, in the IPSL model the temperature anomaly 
spreads over a large area in the Northern Hemisphere, while in the MEP 
model, it concentrates over the area where the ice-sheets are. 

Table [l] compares the global mean surface temperatures obtained using 
the different models, for both Pre-Industrial and Last Glacial Maximum 
conditions. Including the interactive Planck spectrum (version vl com- 
pared to version vO) leads to a slight cooling (0.6 ° C) and a smaller albedo 
sensitivity, while turning off the ozone (version v2 compared to version vl) 
yields a very small warming (0.2 C) and does not change the sensitivity. 
Figure [6] shows the dependence of the global mean surface temperature 
on relative humidity. For simplicity, a horizontally homogeneous relative 
humidity distribution is used. The global mean surface temperature spans 
a wide interval, between approximately 14 C and 24 C. In particular, 
it encompasses the global mean surface temperature obtained with other 
versions of the MEP model and with the IPSL model. 

The latitudinal dependence of surface temperature distributions ob- 
tained from the different modelQ is shown in Fig. [7) for both pre- 
industrial and LGM conditions. When the water vapour feedback is active 
(version v3), the surface temperature is much lower in the polar regions 
than with other versions of the MEP model. For the same reason, the 
response to the albedo change at the LGM is also stronger (Fig. [7] right). 
Globally, the temperature response is approximately 1 ° C stronger than 
in the absence of the water vapour feedback (Table III. 



1 The uniform relative humidity in version 3 is chosen as the mean relative humidity in the 
MEP vO case 
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Figure 6: Solid blue curve: Global mean surface temperature T g as a function 
of relative humidity (with a homogeneous distribution). The horizontal lines 
indicate the temperature obtained by fixing the absolute humidity in the MEP 
model, versions vO (dashed blue), vl (dotted red) and v2 (dashed-dotted yellow), 
and for the IPSL model (green solid line). 





Latitude Latitude 



Figure 7: Left: Surface temperature T g for pre-industrial conditions, for the 
different versions of the MEP model: version vO (solid blue), v2 (dotted red), 
v3 (dashed yellow) and for the IPSL model (solid green). Right: Surface tem- 
perature difference between the Last Glacial Maximum and pre-industrial. 
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Figure 8: Total material entropy production a (solid blue) and global mean 
surface temperature (T g ) (dashed red) as functions of the resolution (number of 
cells). The aspect ratio is maintained equal to 3/4. 

5 The importance of spatial resolution 

In the MEP procedure, it is traditionally argued that maximizing the 
entropy production constitutes a way to represent the effect of small, 
unresolved scales, on the large, resolved scales. In the case of meridional 
heat transport in (dry) planetary atmospheres, the energy is carried partly 
by the mean flow and partly by turbulent fluctuations. Nevertheless, even 
a model accounting for no dynamics at all like the MEP model shown here 
presents reasonable transport curves. For the sake of the comparison with 
the IPSL model, we started with an identical resolution for the GCM and 
the MEP model (Ni at = 72 and Ni on = 96, corresponding to a 3.7 ° x 
2.5° grid). In the MEP model, the resolution is somewhat arbitrary as 
the computational cost is negligible. In the light of the interpretation of 
MEP as a parameterization of small-scale processes, one may naturally 
ask how the results of the MEP model depend on the resolution. 

Figure [8] shows the curves of total material entropy production and 
globally averaged surface temperature obtained with the MEP vO model 
with different resolutions. We keep a constant aspect ratio Ni a t/Ni on = 
3/4 and vary the total number of boxes. Both curves are monotonically 
increasing with resolution. Although there is no explicit representation 
of the dynamics here, the dependence on resolution is very similar to the 
findings of [23] for a GCM. In particular, it shows that the results of the 
MEP model converge when the resolution increase. 



6 Future challenges for MEP climate mod- 
elling 



In this chapter, we have presented a detailed account of how the MEP con- 
jecture can be applied to climate modelling. We have shown how a MEP 
model without ad-hoc hypotheses could be built and we have compared 
its performances in simulating both the pre-industrial and Last Glacial 
Maximum climates with a coupled atmosphere-ocean GCM. The results 
appear to be robust with respect to minor modifications (versions v0-v2) 
of the model. To go beyond these results, we argue that it is necessary 
to account for some feedbacks, and show how to treat them in the MEP 
framework. We stress the importance of the water vapour feedback (ver- 
sion v3) on the surface temperature. Going further would now require the 
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ability to include a water-cycle model in our MEP model. From there one 
may hope to be able to represent clouds in a more robust way than in the 
original Paltridge model. To become a realistic climate model, the MEP 
model would still require important features, like a seasonal cycle, a rep- 
resentation of atmospheric dynamics, a more accurate description of the 
vertical structure, etc, but there are reasons to believe that this would not 
be completely out of reach. This key challenge would have to be taken up 
without sacrificing the original strengths of the MEP model (absence of 
empirical parameterizations and ad-hoc coefficients, rapidity, conceptual 
simplicity). Another major point which would deserve clarification is the 
theoretical basis of the MEP principle (see 7 ). In particular, it would be 
desirable to establish which entropy production should be maximized: is 
it always the material entropy production ? (see for instance |37]). 

If this program could be achieved, the climate modelling community 
would acquire a valuable new tool, in addition to the existing hierarchy of 
models, to improve our understanding of past, present and future climates, 
on Earth and beyond. 
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